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Abstract 

A new variant of the Wolff cluster algorithm is proposed for simulating systems with competing interactions. This 
method is used in a high-precision study of the Lifshitz point of the 3-D ANNNI model. At the Lifshitz point, 
several critical exponents are found and the anisotropic scaling of the correlators is verified. The functional form 
of the two-point correlators is shown to be consistent with the predictions of generalized conformal invariance. 
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Competing interactions are encountered in a 
large variety of physical systems such as, among 
others, magnets, alloys or ferroelectrics. These in- 
teractions may lead to rich phase diagrams with a 
multitude of phases as well as to special multicrit- 
ical points called Lifshitz points (LP). At an LP, 
a disordered, a uniformly ordered and a periodi- 
cally ordered phase become indistinguishable [1]. 
A large number of systems (magnets, ferroelec- 
tric liquid crystals, uniaxial ferroelectrics, block 
copolymers) have been shown to posses an LP. 

This kind of systems may be mimicked by spin 
models with competing interactions. The simplest 
of these models is the well-known ANNNI (ax- 
ial next nearest neighbour Ising) model [2,3]. Be- 
cause of the presence of the competing interactions. 
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previous Monte Carlo studies of the transition(s) 
from the disordered high-temperature phase to the 
ordered low-temperature phases exclusively used 
single-spin flip algorithms [4,5]. Due to the criti- 
cal slowing-down inherent to these methods, those 
numerical studies were limited to rather small sys- 
tem sizes. However, large systems are needed for 
the precise computation of critical quantities, such 
as critical exponents or critical correlators. 

In this contribution, we present a new variant 
of the Wolff cluster algorithm [6] specifically de- 
signed for the simulation of systems with compet- 
ing interactions near criticality. Similarly, we also 
generalize the recently proposed method of Evertz 
and von der Linden for the computation of corre- 
lation functions of infinite systems [7] to systems 
with competing interactions. As an example, we 
present results obtained at the LP of the three- 
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dimensional ANNNI model, and we shall concen- 
trate on the LP critical exponents and critical cor- 
relators. In particular, our data for the scaling of 
the two-point correlation functions are consistent 
with the theoretical predictions of generalized con- 
formal invariance [8,9]. 

Problems coming from critical slowing-down en- 
countered when using local Monte-Carlo dynam- 
ics are alleviated by using non-local methods, such 
as the Wolff cluster algorithm [6]. For the Ising 
model with only a nearest-neighbour coupling J, 
this algorithm may be described as follows: one 
chooses randomly a lattice site, the seed, and then 
builds up iteratively a cluster. If i is a cluster site 
(with spin Si) a site j neighbouring the cluster site 
i will be included into the cluster with probabil- 
ity p = i (1 + sign (s,sj)) (1 - exp [-2 J/(fcsT)]). 
One ends up with a cluster of spins having all the 
same sign, see Fig. la, and which is flipped as a 
whole. This kind of same-sign clusters is obviously 
not adapted to our problem because of the com- 
peting interactions. 

The modified cluster algorithm we propose for 
simulating spin models with competing interac- 
tions is in the following presented for the ANNNI 
model. A generalization to other models is straight- 
forward. The Hamiltonian of the 3D ANNNI model 
is, with Sxyz = ±1, 

H = — J Sxyz {s(x+l)yz + Sx{y+l)z + ^2:1/(2+1)) 
xyz 

+ K J SxyzSxy(z+2) (1) 
xyz 

whereas J > and k > are coupling constants. 
In the axial z-direction competition between ferro- 
magnetic nearest-neighbour and antiferromagnetic 
next-nearest-neighbour couplings occurs, leading 
to a rich phase diagram [10]. 

Our proposed modified cluster algorithm starts 
with a randomly chosen seed and builds up iter- 
atively a cluster. Consider a newly added cluster 
lattice site i with spin Si . A lattice site j with spin 
Sj nearest neighbour to i is included with proba- 
bility pn = p, whereas an axial next-nearest neigh- 
bour site k with spin Sk is included with probability 
Pa = i(l-sign(s,Sfe))(l-exp[-2jK/(fcsr)]). 
Thus the final cluster, which will be flipped as a 
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Fig. 1. Typical clusters (gray spins) obtained (a) by the 
Wolff algorithm and (b) by our modified algorithm for 
systems with competing interactions, here shown in two 
dimensions for simplicity. The competition takes place in 
the axial direction. The spin enclosed by a circle is the seed. 

whole, contains spins of both signs, as shown in Fig. 
lb. At variance with the traditional Wolff method, 
the flipped spins are not necessarily connected by 
nearest-neighbour couplings. Ergodicity and de- 
tailed balance are proven as usual. This algorithm 
works extremely well in the paramagnetic phase, 
in the ferromagnetic phase and in the vicinity of 
the LP. It has not yet been subjected to stringent 
tests in the modulated region where large free en- 
ergy barriers make simulations very difficult (we 
point out recent progress achieved in microcanoni- 
cal simulations of finite ANNNI systems [11]). Note 
that our algorithm differs from the same-sign clus- 
ter algorithm proposed by Luijten and Blote [12] 
for the simulation of systems with long-range in- 
teractions. 

In an interesting work, Evertz and von der Lin- 
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den [7] have proposed a new numerical scheme, 
based on the Wolff cluster algorithm, for the direct 
computation of two-point functions of an infinite 
lattice. In their approach, each cluster is started 
at the same lattice site and not from a randomly 
chosen site. After a certain number of cluster flips, 
the two-point functions can be measured within a 
certain region, this region getting larger with in- 
creasing number of flips. With this method, corre- 
lation functions of the infinite lattice are obtained 
at temperatures above the critical temperature [7] . 
At the critical temperature, finite-size effects are 
greatly reduced as compared to more traditional 
approaches. We have generalized the method of 
[7] to systems with competing interactions by re- 
placing the usual Wolff clusters with our modified 
clusters. This enables us to compute with unprece- 
dented precision two-point correlators in spin sys- 
tems with competing interactions. 

In the following, we report results of a large- 
scale Monte Carlo study of the uniaxial LP en- 
countered in the 3D ANNNI model where the new 
algorithms introduced above have been used with 
great success [9] . The whole study took the equiva- 
lent of more than two CPU-years on a DEC alpha 
workstation. For the investigation of LP proper- 
ties a precise location of this point is mandatory. 
One possibility to locate this point is given by the 
structure factor S{q) = ^{Sq Sr) exp{iqz) as the 

r 

transition between uniformly ordered and period- 
ically ordered phases shows up as a shift of the 
maximum of S{q) from q = to a non-zero value 
[4]. Having found a reliable value for kl, the LP 
critical temperature is then obtained by standard 
methods. In Table 1, we compare the location of 
the LP as obtained from our data with previous 
estimations. Note that Kaski and Selke [5] did not 
attempt an independent determination of kl but 
merely determined Tl for a k in close vicinity to 
Kl as obtained from a high-temperature (HT) se- 
ries expansion [13]. Prom Table 1 the increase in 
precision is evident. 

The uniaxial LP is a strong anisotropic critical 
point where the correlation lengths parallel and 
perpendicular to the axial direction diverge with 
different critical exponents: ~ |T — T^j"''!!'^. 
The value of the anisotropy exponent 9 = iy\\/i'± is 



Tabic 1 

Estimated location of the LP in the 3D ANNNI model de- 
fined on a cubic lattice, as obtained from high-temperature 
expansions (HT) and Monte Carlo simulations (MC). 







Tl 


HT [13] 


0.270 ± 0.005 


3.73 ±0.03 


MC [5] 


0.265 


3.77 ±0.02 


present work 


0.270 ± 0.004 


3.7475 ± 0.005 



Table 2 

Critical exponents at the LP of the 3D ANNNI model, as 
obtained from Monte Carlo simulations (MC) and renor- 
malized field theory (FT). The numbers in brackets give 
the estimated error in the last digit (s). 





a 




7 


(2 - a)/7 


/3/7 


MC [5] 




0.19(2) 


1.40(6) 




0.14(2) 


FT [14] 








1.27 


0.134 


FT [17] 


0.160 


0.220 


1.399 


1.315 


0.157 


present work 


0.18(2) 


0.238(5) 


1.36(3) 


1.34(5) 


0.175(8) 



in good approximation equal to 1/2 [14]. In order 
to take into account the special finite-size effects 
coming from the anisotropic scaling at the LP [15], 
large systems of anisotropic shape with L x L x N 
spins with 20 < L < 240 and 10 < < 100 have 
been simulated. 

Our estimates for the LP exponents a, P, 7 are 
given in Table 2. For the first time, these ANNNI 
model exponents are computed independently. 
They were found by investigating effective expo- 
nents which yield the critical exponents in the 
limit T — !■ Tl [9,16] (provided that finite-size 
effects can be neglected). Our error bars take into 
account both the sample averaging and the un- 
certainty in the location of the LP. The reliability 
of our data is illustrated by the agreement of the 
independently estimated exponents a, (3 and 7 
with the scaling relation a + 2/3 + 7 = 2 up to 
« 0.8%. Furthermore, the agreement with a recent 
two-loop calculation [17] is remarkable. 

Fig. 2 summarizes our data obtained for the LP 
spin-spin correlator C(r_L,r||) = (5'rj_,r|| ■S'o.l.o)- It 
shows the scaling function $(zt) related to the cor- 
relator by 

C(r^,r||) = r-^^/%(ri/r||) (2) 
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Fig. 2. Scaling funtion versus v}-!"^ = (y/r±/r^^)^^^ 

for K = 0.270 and T = 3.7475. The different symbols cor- 
respond to different values of rj^. Inset: comparison of the 
full data set of 1.7 X 10'* points for the scaling function 
<!>(«) (gray points) with the analytical prediction resulting 
from reference [8] (full curve). 

and the scaling dimension x ~ (2 + 0)/(2 + 7//?) 
[8,9]. In this analysis, we use 6 — 1/2, since the 
small deviations from this value obtained in recent 
field-theoretical calculations [14,17] are not yet dis- 
tinguishable from the purely numerical errors in 
our data. The data shown in Fig. 2 have been ob- 
tained after more than five million cluster updates 
for a system with 200 x 200 x 100 spins. They per- 
mit a nice visual test of the data collapse and es- 
tablish scaling, for the first time also at the LP. 

The inset of Fig. 2 is a direct comparison of our 
numerical data with the theoretical prediction re- 
sulting from a generalization of conformal invari- 
ance to strong anisotropic criticality [8] . Applied to 
the LP of the 3D ANNNI model and assuming 9 = 
1/2, the scaling function i}{v) = u~^^/^$(l/w) 
should satisfy the differential equation 

where ai is a constant. It turns out [8,9] that the 
functional form of only depends on a single 
universal parameter p which can be determined 
from the Monte Carlo data by analyzing ratios of 
moments of $. As shown in [9], different moment 
ratios yield consistent values for this parameter. 
The mean value of p obtained in this way has been 



used in the inset of Fig. 2. The nice agreement be- 
tween theory and numerical data provides strong 
evidence for the applicability of the generalized 
conformal invariance to the strongly anisotropic 
critical behaviour realized at LPs. A similar agree- 
ment between numerics and theory has been ob- 
tained for the energy-energy correlation function. 

In conclusion, we have presented a new variant of 
the Wolff cluster algorithm designed for the study 
of systems with competing interactions. Applied to 
the LP of the 3D ANNNI model, the anisotropic 
scaling of the two-point correlators could be con- 
firmed directly for the first time. The form of the 
scaling function is in agreement with the hypothe- 
sis of generalized conformal invariance. 
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